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We present a phenomena-logical analysis of the cos 20 asymmetry recently measured by the COM- 
PASS and HERMES collaborations in unpolarized semi-inclusive deep inelastic scattering. In the 
kinematical regimes explored by these experiments the asymmetry arises from transverse-spin and 
intrinsic transverse-momentum effects. We consider the leading-twist contribution, related to the so- 
called Boer-Mulders transverse-polarization distribution hi (x, fcf.), and the twist-4 Cahn contribu- 
tion, involving unpolarized transverse-momentum distribution functions. We show that a reasonably 
CN) , good fit of the preliminary data sets from COMPASS and HERMES is achieved with a Boer-Mulders 

function consistent with the main theoretical expectations. Our conclusion is that the COMPASS 
and HERMES measurements represent the first experimental evidence of the Boer-Mulders effect 
in SIDIS. 
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I. INTRODUCTION 



43 

Among the various obscrvablcs related to the transverse momentum and the transverse spin of quarks (for reviews, 
Ph| see 0,0), the azimuthal asymmetries in unpolarized semi-inclusive deep inelastic scattering (SIDIS) at moderate Pt 
have recently attracted a large experimental and theoretical attention (see, for instance, [3l4r2j]). In particular, the 
cos 20 asymmetry potentially represents a fundamental source of information on the Boer-Mulders function h^(x, kr) 
[l3j , which measures the transverse polarization asymmetry of quarks inside an unpolarized hadron. 

In a previous paper which predated the experimental measurements, we presented a systematic study of the 
cos 20 asymmetry in unpolarized SIDIS, taking into account both non perturbativc and pcrturbative contributions. 
We found that (cos 20) is of the order of few percent, but our main prediction was that, due to the Boer-Mulders 



effect, the asymmetry in ir~ production should be larger than in ir + production. The recently released COMPASS 
■ [HI, [HI and HERMES data [13] confirm this expectation (the cos 20 asymmetry was previously investigated at high 
\ Q 2 , where it is dominated by perturbative QCD effects [13 ]. by the EMC [HI and ZEUS [l9[ experiments). 

The purpose of the present work is to perform a phenomenological analysis of the recent (cos 20) measurements 
' in order to check their mutual compatibility and to extract some information about the Boer-Mulders function. The 
available data do not allow yet a full fit of h± , with its x and fc^r dependence: so we will simply relate to its chiral- 
even partner, the Sivers function f^j, (describing the distribution of unpolarized quarks in a transversely polarized 
hadron [20l, l2l| ) and limit ourselves to fitting the proportionality coefficient between the two distributions. 
1 We will see that a reasonably good description of both HERMES and COMPASS preliminary data is achieved with 



a Boer-Mulders function close to that used in Ref. 1J] and consistent with various theoretical expectations (impact 



parameter approach [22j , lattice results [23| , large- N c predictions 24 1 and model calculations @, floWT^ . We found 



however that the quality of the fit depends on the assumptions about the average transverse momenta of quarks for 
each measurement. One problem emerges also quite clearly: while the x and z dependencies of the two sets of data 
are satisfactorily reproduced, the Pt behavior of the COMPASS data appears to be incompatible with the HERMES 
corresponding behavior, and hard to understand phenomcnologically. 

II. THE cos 20 ASYMMETRY IN UNPOLARIZED SIDIS 

The process we will consider is unpolarized SIDIS: 

i(£) + P (P) ->■ i'(e') + h(p h ) + x(p x ) . (i) 

The SIDIS cross section is expressed in terms of the invariants 
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FIG. 1: Lepton and hadron planes in semi-inclusive deep inelastic scattering. 



where q = I — I' and Q 2 = —q 2 . The reference frame we adopt is such that the virtual photon and the target proton 
are collinear and directed along the z axis, with the photon moving in the positive z direction (Fig.[TJ. We denote by 
k-r the transverse momentum of the quark inside the proton, by Pt the transverse momentum of the hadron h, and 
by pt the transverse momentum of h with respect to the direction of the fragmenting quark. All azimuthal angles 
are referred to the lepton scattering plane (we call <j> the azimuthal angle of the hadron h, see Fig. [TJ. We follow the 
conventions of Ref . [25( . 

We are interested in the low- Pt region, Pt < 1 GcV. In Ref. [l4| we showed that the perturbative contribution 
to (cos2</>) [26H29| is negligible in this domain. Therefore, in the present analysis we only take into account the 
non-perturbative effects, related to the transverse momentum of quarks. For small transverse momenta it has been 
formally proven that SIDIS can be described in terms of Transverse Momentum Dependent (TMD) distribution and 
fragmentation functions [3(| ■ We will work at tree level, that is in the parton model generalized to include transverse 
momenta of quarks. At leading twist, there is just one contribution, involving the Boer-Mulders distribution hi 
coupled to the Collins fragmentation function [3l[, which describes the fragmentation of transversely polarized 
quarks into polarized hadrons. Concerning higher-twist terms, it is known that there is no twist-3 contribution to 
(cos 2 <b) 13211 . whereas at twist-4 (i.e., at order k T /Q 2 ) the asymmetry receives a contribution from the so-called Cahn 
term [33l. l34j. arising from the non-collinear kinematics of quarks. 

The </>-indcpendcnt part of the SIDIS differential cross section reads 



where ff(x, k T ) is the unintegrated number density of quarks of flavor a and Di(z,p T ) is the transverse- momentum 
dependent fragmentation function of quark a into the final hadron. 

At leading-twist the only fey-dependent term contributing to the cos 2<f> asymmetry contains the Boer- Mulders 
distribution hi coupled to the Collins fragmentation function of the produced hadron. This contribution to the 




(3) 



cross section is given by [l3j 



where Mm is the mass of the nucleon, Mh is the mass of the produced hadron and h = Pt/Pt- 
The twist-4 Cahn contribution has the form 




(4) 
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x J d 2 k T J d 2 p T S 2 (Pt - zk T - pr) 



x 2(kT '^ 2 4 ft(x, D?(*,P§0 cos 20. (5) 

Notice that the Cahn term is only a part of the total, still unknown, twist-4 contribution. 
The asymmetry determined experimentally is defined as 

2 (cos 20) =2^^. (6) 

The integrations are performed over the measured ranges of x, y, z, with a lower cutoff P™ ln on Pr, which represents 
the minimum value of Pr of the detected charged particles. 

Using the expressions above, the numerator and the denominator of © are given by 



j dacos2 ^=^d //// T l ^-y){^\fiM + lB[kt > ,H^]} ) 

TV 

Q 



(7) 

H~ J J J J * 
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d(7= 2™L 1 I I I I ^el X [l + (l-yf]C[f?,D<t], (8) 



where 



P™ ax px 2 py 2 PZ2 

AP T Pr / dx dy dz (9) 

Jxi Jyi Jzi 



and (x is the angle between Pr and kr) 

A[f?, Dl\ = J d 2 k T j d 2 p T S 2 (P T - zk T - p T ) 



2(k T .h) 2 k T cqs20 

Q 

' , [ 2 \ 2fc 2 cos 2 X -fc 2 
dfc T fc T / dx 



xf?(x,k T )D<t(z,\VT-zk T f), (10) 



S^ , i^ ] = J d 2 k T /" d 2 p T c5 2 (P T - zk T - p T ) 



zM N M h 

a, u 'a k T + {P T /z)k T cosx-2/c 2 cos 2 X 



Jo MM h 
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x fc£) ff^, |P T - zk T | 2 ) , (11) 

C[ft,D1] = J d 2 k T y'd 2 p T ^(P T -zk T -p T )/ 1 a ( a; ,fc 2 )^(z,p 2 1 ) 

/•OO /•27T 

= / dk T k T / d X A a (x,fc 2 )^(z,|P T -zk T | 2 ). 
Jo Jo 



(12) 



III. PARAMETRIZATIONS OF DISTRIBUTION AND FRAGMENTATION FUNCTIONS 

To calculate the azimuthal asymmetries we need first of all the /cr-dependent unpolarized distribution functions, 
which we assume to have a Gaussian behavior in kr'- 

p — k T / (k T ) 

f!(x,k T )=m X ) _ /u2 . . (is) 
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= -0.35 


A d = 


0.90 
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0.24 




= -0.04 


A 3 = 


0.40 


As = 


-1 


Otu 


= 0.73 


ad = 


1.08 


a S ea 


= 0.79 


P 


= 3.46 


Ml = 


= 0.34 (GeV/c) 2 







TABLE I: Parameters of the Sivers function used in Eqs. (|16I17|I 



The Gaussian dependence of the transverse-momentum distribution functions is supported by a recent lattice study 
[35l ] . The integrated unpolarized distribution functions ff are taken from the GRV98 fit (36| • 

The available data on (cos 20) do not allow a full extraction of the Boer-Mulders function. Thus we simply take 
hj~ to be proportional to the Sivers function f^,, 

h± 9 (x,k*) = \ <l f^(x,k 2 T ), (14) 

with a coefficient X q to be fitted to the data. Various theoretical arguments (based on the impact-parameter picture 
[22I ] , on large- N c arguments [HI , and on model calculations [Tol - [l^ | ) suggest that that the u and d components of hj~ , 
at variance with f^p, should have the same sign and in particular be both negative (which means that Ad should be 



negative). This is indeed what we find in our analysis. Moreover, the impact-parameter approach 22( combined with 
lattice results (23| predicts a u component of larger in magnitude than the corresponding component of and 
the d components of hi and f^ T with approximately the same magnitude (and opposite sign). 

We parametrize the Boer-Mulders function using the Ansatz (fT4"]) and taking the Sivers function from a fit to 
single-spin asymmetry data [37| . Thus we set 



h 1 q (x,kT) = A g / 1 J J ?(x,fc|) = X q p q (x)r](kT) fi{x,Hj>) > ( 15 ) 



where 



p q (x) = A q (1 - xp K + M^ f (16) 

a q q b q q 

r,(k T ) = VTe^e-^. (17) 

Here Mp is the proton mass, A q , a q , b q and Mi are parameters determined in [37} (see Table H]). Being a quark spin 
asymmetry, f^, must satisfy a positivity bound, which is automatically fulfilled by the parametrization of Ref. [37J. 
Notice that the Sivers function parametrization, as defined in Ref. [33], is: A N f q (x, k±) = —2-^ f^(x 7 k±). 

Concerning the antiquark Boer-Mulders distributions, the SIDIS (at least, the present ones) are not able to constrain 
them. Thus we simply take the Boer-Mulders antiquark distributions to be equal in magnitude to the corresponding 
Sivers distributions and both negative. Note that the Drell-Yan measurements of the cos 2(f> asymmetry [H, [33 would 
in principle give information about the antiquark sector [40l l4l| , but most of the present data seem to be explainable 
in terms of perturbative QCD [42], [43[ . 

Let us now turn to the fragmentation functions. We distinguish their favored and unfavored components, according 
to the following general relations 

Ar+/u = Ar+/d = Ar-/d = Ar-/tt = Aav (18) 
Dir+/d = Dir+fu = Dir-fu = Ar-/d = D n ±/ S = D n ±/ S = D un f, (19) 

The pT-dependent unpolarized fragmentation D±(z,pj,) is assumed to have the form 

p -Pt/<Pt> 

D 1 (z,p 2 T ) = D^z) , (20) 

again with a Gaussian behavior in p^- The integrated fragmentation function Di(z) is taken from the the DSS fit 
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TABLE II: Parameters of the favored and unfavored Collins fragmentation functions [45 



For the Collins function we use the parametrization of [451 ]. based on a combined analysis of SIDIS and e + e data: 



with 



p c q {z)=A c qZ i{l-z) 5K 



rf{jp T ) 



M c 



(21) 

(22) 
(23) 



We let the coefficients A q to be flavor dependent (q = u,d), while all the exponents 7,(5 and the dimensional 
parameter Mc are taken to be flavor independent. The parametrization is devised in such a way that the Collins 
function satisfies the positivity bound (remember that is essentially a transverse momentum asymmetry). The 
values of the parameters as determined in the fit of Ref. [45( are listed in Table [TTJ 

The remaining crucial ingredients to be considered are the average values of k T and p T . Notice that the following 
kinematical relation holds between the transverse momentum of the produced hadron and the transverse momenta of 
quarks: 



(P 2 t) 



(p 2 T ) + z 2 (k 2 T ) 



(24) 



Due to the limitations of the present data sets it is not possible to treat (k T ) and (pj.) as two additional parameters 
to be determined by the fit. We have to make some assumptions about them. 



In our main fit (hereafter called Fit 1) we take (k 2 
unpolarizcd SIDIS cross section performed in Ref. 



and (jp T ) from the analysis of the azimuthal dependence of the 



(k T ) = 0.25 GeV 2 



(Pt) 



0.20 GeV 2 . (25) 

given by their own analysis of 
with (k 2 ) = 0.18 GeV 2 and a 



We also tried another fit ("Fit 2"), using for HERMES the values of (k T ) and {p%_ 
the Pt spectrum, which turns out to be reproduced by Monte Carlo calculations [4t 
^-dependent transverse momentum of the fragmenting quark, {jp T ) — 0.42z 37 (l — z) 54 GeV 2 . In the z range of 
interest this is very well approximated by (p 2 ^) ~ 0.20 GeV 2 . Thus in our Fit 2 we choose for HERMES: 

(k 2 ,) = 0.18 GeV 2 , (p 2 T ) = 0.20 GeV 2 . (26) 

We have no similar information for the COMPASS measurement and for their data we still use in Fit 2 the values 
([25"]) . Therefore Fit 2 is characterized by a (k T ) which is different for the two sets of data. As we shall see in the next 
Section, Fit 2 turns out to be significantly better than Fit 1. 



Finally, concerning a possible flavor-dependence of (k T ) 



47] . we showed in our previous paper 14| that it hardly 



affects the results, hence we shall not take it into account here (it should also be remarked that the experimental 
evidence for a flavor-dependent (k T ) is still far from being established). 
In summary the assumptions of our fits are: 



h 1 q (x, kj.) = X n /jjQe, k 2 ) for 
functions as given in Ref. [33]. 



and d quarks while h^ q (x,k^) = — \f 1 jl(x ) k T )\ for sea quarks, with 



• The Collins functions H 1 9 (z,pt) is taken as in Ref. (451 ] . 

• Gaussian transverse momentum distribution is assumed for all the distribution/fragmentation functions. 
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FIG. 2: Our Fit 1 to the HERMES preliminary proton-target data The dot-dashed line is the Cahn contribution, the 

dotted line is the Boer-Mulder contribution, the continuous line is the resulting asymmetry taking both contributions into 
account. 



• In particular the average transverse momenta for the unpolarizcd distribution and fragmentation functions are 
respectively: 

Fit 1: (k T ) = 0.25 GcV 2 , (p T ) = 0.20 GeV 2 for both HERMES and COMPASS data. 
Fit 2: (k T ) = 0.25 GcV 2 , (p T ) = 0.20 GeV 2 for COMPASS data; 
(k T ) = 0.18 GcV 2 , (p T ) = 0.20 GcV 2 for HERMES data. 



IV. ANALYSIS OF THE (cos 24>) DATA 



Data on the cos 2<f> asymmetry in unpolarized SIDIS at small Pt have been recently presented in a preliminary form 
by both the HERMES Collaboration and by COMPASS [H, EH (while this analysis was in progress the CLAS 
Collaboration at JLab has released some results on (cos2</>), but their conclusion is that the precision of the data 
does not allow obtaining information about the Boer-Mulders function ;48j). The first qualitative evidence coming 
from both COMPASS and HERMES measurements is a larger asymmetry for tt~ production compared to 7r + . This 
difference was predicted in Ref. [l4| to be a signature of the Boer-Mulders effect, which has opposite signs for ir + and 
7r~ (whereas the Cahn contribution, is the same for tt + and 7T~). 

Fitting the HERMES and COMPASS data as explained in the previous Section we find in Fit 1 the following values 
for the coefficients A u and A^: 

A u = 2.0±0.1, A d = -1.111 ±0.001 (Fit 1) (27) 

This implies that hi and hi are both negative. The \ 2 P er degree of freedom of Fit 1 is x 2 /d.o.f. = 3.73. The 
value of Xd corresponds to the saturation of the positivity bound of . Errors on the parameters were calculated 
with Ax 2 = 1. 

Notice that we have excluded from the fit the COMPASS data in Pt , which - as wc will sec below - are clearly 
incompatible with the HERMES Pt data and have a counter-intuitive behavior (it is in fact difficult to envisage a 
transverse- momentum dependence of distribution and fragmentation functions able to describe them). 

The first moments (in k T ) of the Boer-Mulders distributions hf- u and hi, 

ht W (x) = Jd^ T ^- 2 hi(x,k T ), (28) 
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FIG. 3: The first moments of hi" and h± d from Fit 1. 



are displayed in Fig. [3] 

In Fig. [2] we plot A cos2 ' t ' = 2 (cos 20) for 7r + and tt~ production at HERMES with a proton target, with the results 
of Fit 1. The asymmetry is shown as a function of one variable at a time, x, z and Pt] the integration over the 
unobserved variables has been performed over the measured ranges of the HERMES experiment, 

Q 2 > 1 GcV 2 , W 2 > 10 GcV 2 , P T > 0.05 GeV (29) 
0.023 < x < 1.0 , 0.2<z<1.0, 0.3 <y< 0.85 
0.2 < x F < 1 . 

Note that the Boer-Mulders contributions to ir + and 7r~ production are opposite in sign. In fact, we have 

(cos 20)5m ~ e 2 hi«(x) H^(z) + e 2 h± d (x) H^(z) , 

(cos 2<f>y BM ~ e 2 u hi u (x) H^(z) + e 2 hi d {x) H^(z) , (30) 

and, as far as 7J r unf (z) ~ -H^ y (z) [H, one gets different signs for the Boer-Mulders effect for positive and negative 
pions. The combination of the Boer-Mulders term with the Cahn term, which is positive and exactly the same for ir + 
and tt~ (if the £>r-dependence of the distributions is flavor blind) gives a resulting asymmetry which is larger for 7r~ 
than for ir + . 

Fig. [4] shows our Fit 1 to A cos2( ^ h at HERMES with a deuteron target. We have neglected nuclear corrections and 
used isospin symmetry to relate the distribution functions of the neutron to those of the proton. 
The experimental cuts of the COMPASS experiment (which runs with a deuteron target) are: 

Q 2 > 1 GcV 2 , W 2 > 25 GcV 2 , (31) 
0.2<z<0.85, 0.1<y<0.9 

Pt > 0.1 GeV . (32) 

In Fig. [5] we show our fit to the COMPASS data. One clearly sees that the Pt dependence of these data is incompatible 
with the HERMES one and hard to understand theoretically. 

Let us now come to Fit 2. In this case, the coefficients X u and are found to be 

A„ = 2.1±0.1, X d = -1.111 ±0.001 (Fit 2), (33) 

and are very close to those of Fit 1 (again, saturates its positivity bound). Thus, the x-dependence of the Boer- 
Mulders functions is essentially the same in the two fits. However, the \ 2 P er degree of freedom of Fit 2 is significantly 
smaller: y^/d.o.f. = 2.41. In Figs. El [7] and [8] wc show the results of Fit 2 for A cos2 ^ compared to the data. 

In both Fit 1 and Fit 2 we used only statistical and, when provided (HERMES), systematic uncertainties of the 
experimental data. Given the quality of present data and the theoretical uncertainties related to the twist four 
contributions, we did not attempt a more sophisticated analysis of the x dependence of Boer-Mulders functions and 
used the simplified assumption of Eq. 1141 

The main difference between the two fits resides in the Cahn term, which is strongly sensitive to the average value 
of k T . The fact that the data prefer the fit with (fcy)HERMES ^ (k T ) compass seems to indicate that the twist-4 
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FIG. 4: Our Fit 1 to the HERMES preliminary deuteron-target data. The line labels are the same as in Fig. [2] 
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FIG. 5: Our Fit 1 to the COMPASS preliminary data (deuteron target) [l5|, EH- The line labels are the same as in Fig. [2] 



contributions are different in the kinematics of the two experiments. This is a further indication of the fact that the 
cos 20 fit is strongly affected by twist-4 contributions, which are not yet fully known. 

Our analysis shows that, as far as the x and z dependencies are concerned, both the HERMES and COMPASS 
data are fairly well described. The resulting Boer-Mulders distributions of quarks have the expected sign fiol [241 [50l] . 
Moreover, looking at eq. (f2~T)) or ([55)1 . one sees that the Boer- Mulders u distribution is larger by a factor 2 compared 
to the u Sivers distribution, whereas the Boer-Mulders and Sivers d distributions have approximately the same 
magnitude. This is in agreement with the predictions of the impact-parameter approach [22| combined with lattice 
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FIG. 6: Our Fit 2 to the HERMES preliminary proton data. The line labels are the same as in Fig. [2] 
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FIG. 7: Our Fit 2 to the HERMES preliminary deuteron target data. The line labels are the same as in Fig. [2] 



results [HI. 



V. CONCLUSIONS AND PERSPECTIVES 



The present unpolarized SIDIS data on azimuthal cos 2<p asymmetries, although still preliminaries, represent 
clear manifestation of the Boer-Mulders effect. However, they are not sufficient to allow a full extraction of hj~. 
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FIG. 8: Our Fit 2 to the COMPASS preliminary data. The line labels are the same as in Fig. [2] 

particular, twist-4 Cahn contribution turns out to be comparable in size to twist-2 Boer-Mulders contribution so 
measurements of the cos 20 asymmetry at different values of Q 2 are desirable. Measurements at proposed Electron 
Ion Collider [5lll52| at higher Q 2 would allow separation of pure twist-2 from higher twist contributions. To minimize 
the influence of Cahn contribution which is supposed to be flavor blind one can measure cos 20 asymmetry of the 
difference 7r~ — tt + . Moreover, the antiquark Boer-Mulders distributions turn out to be essentially unconstrained. 
Therefore, it would be important not only to have more SIDIS data (the JLab experiments will play in the next 
future a very important role in this respect), but also to explore other processes, like Drell-Yan (DY) production. The 
kinematics of the recent E866/NuSea DY data [H, [39| is such that they are partly dominated by perturbative effects. 
An attempt to extract from these data has been made in Refs. [13, Hi]]. The first moments of quark Boer-Mulders 
functions found in these analysis agree with our results in the relative sign of u and d quark Boer-Mulders functions. 
The magnitudes cannot be easily compared. First of all, as we said, the extraction of hi from E866/NuSea data is 
affected by perturbative effects which in Refs. [4(| Slj, have not be taken into account even at large qx- Second in 
pp and pD DY processes it is possible to extract only the products of quark- antiquark Boer-Mulders functions. In 
order to separate them, in Refs. [13, S3 the positivity bound has been employed, obtaining an allowed range for each 
distribution. However one should add to this range the statistical errors of the fit, which results in a very large final 
uncertainty on the magnitude of h j~ . 

The pp or pD cos 2<j) DY asymmetry, which involves the sea distributions, is very small. A larger asymmetry is 
predicted for pp DY production [53l - l57j , a process to be studied in the next years at the GSI High-Energy Storage Ring 
[58h60I |. Another reaction probing valence distributions is ttN DY, under investigation by the COMPASS collaboration 
in their hadronic program {6lj . 

A combined analysis of all these data will represent a decisive step towards a better knowledge of the Boer-Mulders 
function and of the phenomena originating from it. 
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